Exact results for the entanglement across defects in critical chains 
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We consider fermionic and bosonic quantum chains where a defect separates two subsys- 
tems and compare the corresponding entanglement spectra. With these, we calculate their 
Renyi entanglement entropies and obtain analytical formulae for the continuously varying 
coefficient of the leading logarithmic term. For the bosonic case we also present numerical 
results. 

I. INTRODUCTION 

In critical quantum chains, the entanglement entropy between a section of length L and the 
remainder varies as In L with a prefactor proportional to the central charge c of the model. For a 
review, see If one modifies the interface, this coefficient has been found to vary continuously 
with the defect strength in free particle systems. The effective number of states in the Schmidt 
decomposition then increases as a power of L with a continuously varying exponent. 

For fermionic systems, the problem was first posed by Levine 2| and then investigated numer- 
ically for XX chains [3| and transverse Ising chains [J]. By mapping the problem to that of a 
two-dimensional Ising model with a defect line, the coefficient could later be obtained analytically 
and perfect agreement with the numerical data was found [5j]. Moreover, the parameter entering 
the analytical expression turned out to be simply the transmission amplitude through the defect. 
This holds also for more complicated defects Q|. In a series of recent papers, calculations were also 



7J9]. On 



the bosonic 



done for continuous fermionic systems, and the same coefficient was found 
side, Sakai and Satoh studied a continuum system with a conformal interface between two different 
critical parts with c = 1 [lOj]. This can also be viewed as a uniform system with a defect. The 
continuously varying coefficient found in this case has a close relation to the fermionic one, but the 
detailed connection has remained unexplored so far. 

The purpose of this note is two-fold. Firstly, we want to show how one can treat the bosonic 
case in complete analogy to the fermionic one. Thus we consider a system of coupled oscillators 
which is the lattice version of the system studied in [ljj]. For this system, we derive an expression 



2 



for the single-particle spectrum in the reduced density matrix p which is formally very similar to 
that in the fermionic case. However, it does not have a gap, since the defect does not break the 
criticality. This clarifies the relation between the two problems at the level of the RDM spectra. 
We check that one recovers the previous result for the von Neumann entropy in this way and also 
compare with numerics. 

Secondly, with the spectra at hand, we show that not only the von Neumann entropy S = 
— tr (plnp), but also the Renyi entropies 

S n = T ^\ntT{p n ) (1) 

can be calculated asymptotically in closed form for all integer n. In both cases, the coefficients K n 
in 

S n = n n lnL (2) 

turn out to be (sums of) elementary functions for n > 1 and in this sense simpler than the von 
Neumann coefficients. Qualitatively, they all vary in a similar way for each of the two cases. 

In the following, we first recapitulate the fermionic results in section 2. Then, in section 3, we 
study the oscillator chain and its RDM spectrum. In section 4 and 5 we present the calculations 
for the Renyi entropies and show the resulting functions. In section 6 we sum up our findings and 
in the appendix we present a derivation of the bosonic spectrum from the transfer matrix of a 
two-dimensional Gaussian model. 

II. SETTING AND FERMIONIC RESULTS 

We consider open chains of length 2L with a defect in the middle and the entanglement between 
left and right halves. The geometry is shown in fig. 1 for a bond defect. 
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FIG. 1: Quantum chain with a bond defect. 



For free-particle systems, the reduced density matrix for a subsystem can be written 

1 u 



(3) 
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where H is again a free-particle Hamiltonian, see TJ|. Its single-particle eigenvalues, called 2w/% in 



the following, contain the basic entanglement information. With them, the von Neumann entropy 
S is given by 

fc fc 

and the Renyi entropies are 

S n = j^— [ ± M 1 ± e" 2nWfc ) T n ^ ln(l ± e" 2aJfc )] (5) 
fc fc 

where the upper (lower) sign refers to fermions (bosons). 

For a transverse Ising model, these eigenvalues were determined in [jj from the excitations in 

the transfer matrix of an Ising strip of width InL. The result was that 

chw fc = - che fc (6) 
s 

Here the are the values without the defect and vary linearly with k for InL 3> 1, while s 
measures the defect strength. For a bond defect, where the coupling is changed from 1 to t, it is 
given by s = sin(2 arctan t) = 2/(t + 1/t). For an XX chain, T = s 2 is the transmission coefficient 
through the defect at the Fermi level. The dispersion relation ([6]), shown in fig. 5 of [5|, describes 
a spectrum with a gap induced by the defect and encodes the entanglement properties for a large 
system. In particular, the entanglement decreases as s becomes smaller. This was investigated in 
[5] for the von Neumann entropy. 

The functional relation ([6]) can also be written in another form. In terms of the quantities with 
and without the defect, 

C' k = l/(e 2 ^ + l), Cfc = V(e 2£fc + 1) (7) 

it becomes 

C fc (i-Cfc) = * 2 Cfc(i-Cfc) (8) 

Such a relation was found recently in the study of quantum wires, i.e. for continuous systems with 



a localized, scale- free scattering potential [a, [9(. In these calculations, one considers the overlap 
matrix A of the occupied single-particle states in the subsystem (l^ j and finds that 

A'(l - A') = s 2 A(l - A) (9) 

for arbitrary particle number N. This relation extends to the eigenvalues a/% which are then used to 
determine the entanglement entropies. However, for free particles the non-trivial are the same 
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as the eigenvalues Ck of the correlation matrix C which give the u>k and £/% via (|7|). The relation for 
O'k in ^ 2] is therefore the same as (JSj) an d the calculations, although they proceed in a different 
way, have actually the same basis and thus lead to the same results. 



III. BOSONIC CHAIN 



A quantum chain realizing the system studied in [10J consists of 2L harmonic oscillators coupled 
by springs, where the spring constants and the masses are different in both halves, but have the 
same ratio, see [15]. The spring in the centre has to be chosen properly. The Hamiltonian is 

2L , , o9 ^ 2I/-1 

x n - x n+ i) 2 (10) 



/ 1 d 2 1 \ 1 

H = ( — — — -\ — m„ fin x 2 ) -\ — D„ (x„ 

^ V 2m n dxl 2 n J 2 ^ V 

n=l x " 7 n=l 



and we set 




, ... , n < L 

D n = m n ={ n (11) 

n > L 

while the central spring is assumed to have 

A rescaling of the coordinates u n = ^Jra n x n then makes H homogeneous up to the springs at sites 
L and L + 1 

^/i 5 2 i \ i lfK 02 K 02 Kp 



n=l v n 7 n^L 



L^2 

(13) 

In this form, one is dealing with a defect problem, and the defect is completely characterized by 
the parameter 9. For large \9\, i.e. if K% and K% are very different, the chain is cut in the middle. 
The eigenfrequency fin of the oscillators is included to avoid a zero mode, but will be taken small. 
The single-particle eigenvalues in the RDM follow from the matrix C = 2X 2P containing the 



position and momentum correlations in the subsystem 
eigenfunctions (f> m (i) of the dynamical matrix, one has 

2L-1 

f2 m 



111 ]. In terms of the eigenvalues f2^ and 



Qj = ^2 ^-A mn 4> m (i)4> n (j) (14) 

m,n=0 



where the reduced overlap matrix 

L 

A m n = ^2<P m ( y l)M l ) (15) 



1=1 
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comes from taking the product of X and P in the subsystem (chosen as the left half-chain) . 
In the homogeneous system 



„o / 77171" \ 

nl + 2(1 -cos— J, m = 0,...,2L-l 



and 



(t>m(i) 



I (i- 1/2) mir 
cos — m^O, 



Then the matrix A is for ra,n/0 



2L 



sin ^ (m — n) sin ^ (m + n) 




2Z 



sin -£r(m — n) sin ^-(m + n) 
and if one or both the indices are zero one has 



A 



mO 



A 



1 sin 
"" a = V22L sin^m' 



00 



1/2 



Note, that A mra vanishes if m — n 7^ is even. Furthermore, A mm = 1/2 and one can write 



— 2 % + Dij 



(16) 



(17) 



(18) 



(19) 



(20) 



where Dij is the piece of (|14p with the sum restricted tora-n odd. 

For the inhomogeneous system, the eigenvalue equations can be satisfied by choosing the same 
spectrum as in the homogeneous case Q' m = VL m and making the ansatz 



(21) 



oc m (j) m {i) 1 < i < L 
(3m4>m(i) L <i<2L 

Inserting this into the two modified equations and requiring orthonormality, one obtains the con- 
ditions 



The solutions are 



"m/3m = K , a m a n + /3 m /3 n = 25„ 



(22) 



a m = 1 ± tanh 9, (3 m = 1 =F tanh ( 



where the upper (lower) signs refer to even (odd) indices m. This yields 

A mn (l — tanh 2 6) m — n odd 



A' 



A mn (l + tanh#) 2 m = n even 
A mn (l — tanh6*) 2 m = n odd 



(23) 



(24) 



6 



Finally, one obtains 

Cj. = -(1 + tanh 2 0)8ij + (1 - tanh 2 Q)D^ 

2 ' (25) 

= (1 - tanh 2 6)Cij + tanh 2 B <% 

which is an exact relation between the matrices with and without the defect and translates to their 
eigenvalues coth 2 uj k and coth 2 e k . Written differently, it takes the form 

shcjfc = - shefc (26) 
s 

where s = l/ch#. This is the analogue of the relation (J6j) and has a striking similarity to it. 
However, while the two relations are identical for large u and e, the lower part of the spectrum is 
different. There is no gap in (j26j) . u approaches zero with slope 1/s as e goes to zero. The defect 
does not make the system non-critical. Full dispersion curves are shown in fig. 2 for several s. 
Also marked are the discrete numerical values which are obtained for a chain of L = 200 sites. We 
also note that with 

n' k = l/(e 2 ^-l), n k = l/(e 2 ^-l) (27) 
the functional equation can be written as 

n' k (l+n' k ) = s 2 n k (l+n k ) (28) 
which is the analogue of (El). Furthermore, one should mention that the overlap matrix (1181) has 

fl Q 

a similar structure as the one in [8j, l9j and also satisfies ([9]). The two calculations are therefore 
closely related. 

As in the fermionic case, the parameter s has the meaning of a transmission amplitude through 
the defect. This can be seen from the 2x2 transfer matrix for the scattering problem. This matrix 
has eigenvalues (e s ,e~ e ) and transmits the odd resp. even functions (|17p by multiplying them with 
these factors. This explains the form (|2ip of the perturbed eigenfunctions. One can also see, that 
only the choice f)12|) leads to a transmission coefficient T = s 2 independent of the wavelength (i.e. 
to a scale-free defect) and to such a simple structure of the problem. 

The relation (|26|) can also be derived by going to two dimensions and studying the transfer 
matrix of a Gaussian model. This is sketched in the Appendix. The relation is also implicit in 
the work of [10], but not really visible. However, one can insert it into the bosonic von Neumann 
entropy (HJ), change to integrals and obtain k = I/ir 2 with 

,= -jfd. mi -•-»*) +jT* ^ m 
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FIG. 2: The relation (|2l)]) for the bosonic single-particle excitations in H for several values of the defect 
strength. The points are the numerical values for chains with L — 200. 



Expressing this in terms of hyperbolic functions, differentiating twice with respect to s and using 
a partial integration gives, with x = she 



I"(s) 



dx- 



x 



arsh(— ) 



o Vx 2 + 1 Vx 2 + s 2 s 

1 



1-s 2 



In s 



(30) 



This is (4.12) in [lOj] and the negative of the corresponding result in the fermionic case A 
further difference appears in the integration over s, since I'(0) = 7r 2 /4 here, while it vanishes for 
fermions. Therefore k(s) contains a term linear in s and reads 



k(s) = —s — kf(s) 



(31) 



where the second part is the fermionic result 

k f (s) = -^2 {K 1 + «) ln (! + s) + (l-s) ln(l - s)] Ins + [(1 + s)Li 2 (-s) + (1 - s)U 2 (s)}} (32) 

with h'12 denoting the dilogarithm. This is, written somewhat differently, the final result in [lol j. 
The function c c g = 6/c(s) is shown in fig. 3 and rises smoothly from zero to 1, while k(s) itself 
varies between and 1/6. 

The figure also contains data points from numerical calculations. In order to see a logarithmic 
behaviour of S in the bosonic case, one has to work in a situation, where the correlation length 
given by 1/CIq is larger than the system size but still finite. In the homogeneous system one can 
then obtain n = 1/6 with high precision (3-4 decimal places) from sizes L = 100 — 500. The 
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FIG. 3: The function c e g = 6k(s) for the bosonic interface problem. The points are numerical results, see 
text. 

calculations are done such that 2LVLq = A is held constant as L is varied. With the defect, the 
dispersion curve shown in fig. [2] becomes rather steep for small s and only a single eigenvalue is 
found on the slope for usual sizes 2L. The results depend sensitively on this eigenvalue and deviate 
from the theoretical curve for small s. We have therefore plotted only data for s > 1/4, obtained 
with A = 0.01, where the agreement is very good. 



IV. FERMIONIC RENYI ENTROPY 



We now turn to the Renyi entropy. In terms of the single-particle eigenvalues, it is given by 
([5]). Converting the sum over k into an integral over e leads to the logarithmic behaviour (J2J) of S n 
with coefficient 

1 1 



7r 2 1 — n 



-K„ 



K n = I n - nh 



where 



roo 

I n = de ln(l + 
J o 



-2nu)\ 



(33) 



(34) 



In terms of hyperbolic functions, the quantity K n is 



K n = de [ln(2chnw) -nln(2cho;)] (35) 
Jo 

To evaluate it, one first writes chnw as a product by using 1.391 of [13j . For even n, this gives 

(36) 



chnuj = Y\ 
k=l 



n/2 r ch 2 W -c? 
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where = sin((2fc — l)-/r/2n) and = cos((2/c — \)ir/2n). One now inserts chu; = che/s and 
takes the derivative with respect to s. This gives, with x = che, 

a ™/ 2 9 2 

'M^H-n-g-*^ ,37) 

The first term is compensated by an identical one from n ln(2ch oo) and one has 

n/2 



d.r 



i Vx^l 



k=l 

2 arcsm(scfcj 



1 1 



Cfc — , (38) 



fc=i ^l-s 2 c| 
which can be integrated to give 

n/2 

-Kn(s) = ~ arcsin 2 (gcfc) (39) 
fc=l 

With k —7- (n/2 + 1 — A;) the Cfc can be changed to and one finally obtains 

n/2 

1 1 v-^ 2/ — 1) 7F n 

K n(s) = > arcsm (ssfc), s fe = sin ), n even (40) 

7r z n — 1 ' 2n 

fc=i 



For odd n, the expression (j36j) is changed into 

(n - 1)/2 r ch^-c| 



ch no; = ch u> Y\ 



s 2 



(41) 



k=l 

but leads to the same result as in (|39p . up to the summation limit. Replacing the again by sine 
functions, one has 

1 1 ( ^ /2 2kn 

K n(s) = -n r V, arcsin 2 (ss fc ), s fc = sin(— — ), n odd (42) 

7r z n — 1 ' 2n 

k=l 

Formulae (|40p and (|42[) give the Av n for all integer n. One sees, that they are all elementary functions, 
namely sums of n/2 or (n — l)/2 arcsin 2 terms. This is in contrast to the limiting cases n — > 1 
(yon Neumann entropy) and n — > oo (largest eigenvalue of p), where the dilogarithm Li2 appears 
5j. The result is particularly simple for n = 2,3 where only a single term is present 

K2 (s) = -s arcsin 2 (s/y/2) (43) 

7T Z 

«3(s) = ^ arcsin 2 (sVE/2) (44) 



The first formula was given before in 
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FIG. 4: Coefficients n n for a transverse Ising chain with a defect. Upper panel: as a function of the 
parameter s. Lower panel: as a function of the defect bond strength t. The n-values are, from top to 
bottom, n = 1/2, 1,2,3,4,8,00. 



Qualitatively, the n n all vary similarly with s, rising from zero for s = (the dissected system) 
to a limiting value 



(45) 



for s = 1 (the homogeneous system) which can be obtained either directly from I n (l) = 7r 2 /24n, 
or by carrying out the simple sums which remain in n n for s = 1. For n = 1, Ki(l) = 1/12 = c/6 
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with c = 1/2. The behaviour near s = is quadratic and given by 

K n(s) = -^^-s 2 , S^O (46) 

47r z n — 1 

As n — > 1, the curvature diverges which signals the s 2 ln(l/s)-behaviour one has in the von Neu- 
mann entropy. 

The n n (s) are shown for several values of n in the upper part of fig. 4. In the lower part, the 
same K n are plotted as functions of the bond variable t. The linear behaviour near s = 1 is then 
turned into a quadratic one near t = 1, and the functions become symmetric under t — > 1/t. 

Fig. 4 includes the results for n = 1 and for n = oo which were found in [5]. Also shown is the 

quantity n n for n = 1/2, which can be obtained from the derivative 

1 dr 1 

K > J S ) = L / ax _L_ 

= ^ ~/f _ s 2 ( ff ~ arcsin g ) ( 47 ) 



which gives 

K i/2( s ) = 7r~9 arcsin s (tt — arcsin s) (48) 

It varies linearly, Ki/ 2 (s) — s/2tt, for small s, has negative curvature everywhere and lies signifi- 
cantly above the n = 1 curve, but otherwise fits into the overall scheme. 

The formulae given above all refer to the transverse Ising model, where c = 1/2. For the XX 
(hopping) model with c = 1, they have to be multiplied by a factor of two. In this case, it is also 
easy to construct a scale-free defect in analogy to the bosonic chain. One only has to supplement 



the modified bond t with site energies ±Vl — t 2 at L and L + 1, respectively. Then Q is satisfied, 
the transmission amplitude is s = t and the K n (s) are the relevant quantities. 

These results can also be obtained for continuum systems by working with the overlap matrices. 
Then the particle number N appears in ([SJ instead of L and the K n are found in the form of an 



infinite scries in die parameter * which can be recognized as thai of the function arcsin" 2 j^,!^]. In 
the cited papers, the notation is somewhat different: 2n n is called and T is used instead of s. 



V. BOSONIC RENYI ENTROPY 

The calculation of S n in the bosonic case is very similar. The quantity K n introduced in (f33|) 
becomes 

[•CO 

K n = - de [ln(2sh noj) — nln(2shw)] (49) 
Jo 
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and the necessary formula for shrew is, for odd re 

(n-l)/2 



sh nw = nshw jQ 
This gives the derivative 



k=l 



1 + 



sh 2 uj 



, Sfc = sin(/i;7r/n) 



(50) 



(n-l)/2 



K'Js) 



fc=i - 70 



dx 



Vx 2 + 1 x 2 + s 2 s 2 



(n-l)/3 

E 

fc=i 



1 — s 2 s 



2 S 2 



and the integrations leads to 

(n-l)/2 



arcsin (wl — s 2 s\) 



(n-l)/2 



(51) 



K n {s) = — arcsin (ssfc)[7r — arcsin (ssfc)] = — [7r 2 /4 — arccos 2 (ss^)] (52) 



fe=l fc=i 
After changing to cosine functions in the arguments, one has 



K n {s) = \—^— - V 
7r re — 1 

fc=l 

In the same way, the case of even re can be treated and leads to 
1 



(n-l)/2 



7T- 



arccos 2 (sc*,) 



c fc = cos( 



(2k - 1)tt, 
2re ' 



re odd 



(53) 



«n(s) 



1 1 
^2 ~1 ] 2 



7T _ 



arccos s 



(n-2)/2 r 2 

fc=i 



arccos (sck) 



c k = cos( 



,2/c71\ 

~2n~' 



re even 



(54) 



These are the analogues of the fermionic formulae (|40p and (|42p and one sees that basically the 
arcsin 2 has been replaced with arccos 2 . The value for s = 1 is now 



1 , In 
«• (1) = —(1 + - 



whereas the behaviour for s — > is linear 

1 1 



K n (s) 



7T 



2vr 2 re-l C ° t( 2^ )s ' 



(55) 



(56) 



The limit n — > 1 can be taken and gives the slope 1/4 found already in (|3ip . Again, the cases 
re = 2, 3 give the simplest formulae 



1 




4 




1 - 




8 




1 




4 




1 - 
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arccos 2 (s) 



arccos 



(57) 
(58) 



In fig. 5, the functions K n (s) are shown for six different values of re. They all rise in a rather 
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s 



FIG. 5: Coefficients K n (s) for the bosonic defect. The n-values are, from top to bottom, n = 1, 2, 3, 4, 8, oo. 

smooth way with negative curvature. Included is also the limit n — > oo, for which one only has to 
calculate the integral I\. This is 

/•oo 

h = - de pn(2shw) - u] (59) 
Jo 

and has the derivative 

I{(s) = ^ln(l + s) (60) 
which upon integration gives a dilogarithm and 

Koo = \U 2 {-S) (61) 

Thus Kqo, which describes the scaling of the largest eigenvalue of p, is a simple, but non-elementary 
function, as in the fermionic case. 

VI. CONCLUSION 

We have considered defects in critical chains of free particles which can be varied in such a 
way that one can go continuously from a homogeneous system to one cut in two pieces. In the 
RDM, this leads to a characteristic rise of the single-particle eigenvalues by which the entanglement 
across the defect becomes smaller. Both for fermions and for bosons, the change of the eigenvalues 
is described by a simple functional relation. This allows to obtain closed expressions for the 
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entanglement entropies in the asymptotic region. In this sense, one is dealing here with a fully 
soluble problem. 

The bosonic chain, modeling a system with a conformal interface, is somewhat subtle. In order 
to avoid a zero-energy mode, which would spoil the correlation function approach and does not 



contribute to the logarithmic term anyway 



10l | . one has to work slightly off-critical. The numerical 



extrapolations for the available sizes then do not give the same perfect agreement with the analytical 
results as for the fermionic defect. Nevertheless, they reproduce them over most of the parameter 
space. 

An interesting point is the validity of the functional relation (|26p also away from criticality. This 
offers the possibility to obtain exact results also in this case, since the £/% are explicitly known and 
equidistant in the infinite system [ill. lit]]. In the critical region, one then finds the same behaviour 
as in ([2]) with L replaced by the correlation length £. One can also see that one needs quite large 
values of £ to observe the exact value of K n (s) if s becomes small. 

Finally, one should mention that simply reducing the central spring constant in the oscillator 
chain does not lead to a varying n n . The eigenvalues ujk then increase in a similar way as here 
and the entanglement entropy becomes smaller, but the asymptotic variation with L remains 
unchanged. 



Acknowledgments 

We thank Pasquale Calabrese for stimulating correspondence on the topic and Jens Eisert for 
an interesting discussion. V.E. acknowledges financial support by the ERC grant QUERG. 

Appendix: Transfer matrix in the Gaussian model 

We want to show here that the relation (1261) can also be obtained from the transfer-matrix 



for the fermionic case. 



excitations in the associated two-dimensional model, as done in 

Consider a Gaussian model on a square lattice with variables cf>, — oo < <j) < oo, and coupling 
constant K such that K (0 — 4>') 2 /2 is the energy of neighbouring sites. The symmetrized row 

1 /2 1 /2 

transfer matrix W = V l V2V1 consists of the contribution from vertical (Vi) and horizontal 
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(V 2 ) bonds. These are given by 



V 1 = exp( l -K* £^f), V * = exp( —K J> n - n+1 ) 2 ) (62) 



n ' " n 



where K* = 1/K. After a Fourier transformation with open boundaries and momenta q, they 
become 

V 1 = exp( l -K* £^), ^ = ex P ( ^) (63) 
q q <? 

with ft 2 . = 2(1 - cos q) = 4sin 2 (g/2). These can be expressed in terms of creation and annihilation 
operators b\,b q for oscillators with mass 1 and frequency fl q . Then 

V l = exp( \k* V q (b q - b\f) , V 2 = exp( —K^ + b\f) (64) 
q q 

1/2 

In order to obtain If as a single exponential, one forms Heisenberg operators with V 1 and Vi- 
These are 



(65) 



b\) \ -2A 2 1 - 2A 2 

with the notation Ai = K*Q q /4: and A 2 = KQ, q /A. This gives for W 

.i5)^-^^(;)-(;;)(ji 

The antisymmetric matrix on the right has determinant 1 and thus eigenvalues e ± " lq where ch j q = 
(a + d)/2 is half the trace. Explicitly, 

ch 7 , = l + ^, or sh ( 7 ,/2) = fl q /2 (67) 

Introducing new boson operators via 

the row transfer matrix then becomes 



q 



(69) 
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FIG. 6: Strip geometry defining the transfer matrix Wu 



The parameter in the transformation (|68p is given by 



exp(-40 a ) = -^(l + S) 



(70) 



and depends on the value of K, in contrast to 7 9 . 

Consider now a strip as shown in fig. 6 where the three sections have coupling constants K±, 
K 2 and K\. The total transfer matrix then is Wtot = W^^W^W^ 2 and the relation analogous 
to (1661) contains the matrix 



M tot = u(d 1 )- 1 E M ' 2 u{e 1 )u{d 2 )- 1 E M u{e 2 )u{6 l )- l E M i 2 u{e 1 ) 



(71) 



where E is diagonal with entries e^ 79 and the index q of the 9 has been suppressed. But the 
products of the LPs are just U{9\ — 82) = U{9) where 



exp(2#) = |i 

is independent of fL. Using cyclic permutation, the trace is that of the matrix 



(72) 



Mtot = U(9)E M U(6)- L E 



(73) 



and calling the eigenvalues e ±2ulq one finds 



ch2wg = ch 2 (M 7(? ) + ch(20) sh 2 (M 7(? ) 



(74) 



or equivalently, writing M~f q = e q and ch = 1/s 



sh.UJ q = - sh.E q , 



(75) 



This is the relation (|26|) . Due to (|72p. the parameter s is defined in the same way as for the chain. 
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Remarks: 

(i) Multiplying the symmetrized transfer matrices W\ and W2 leads to particular vertical bonds 
at the interface: Kq = (K{ +K%)/2, i.e. 1/K = (1/Ki + l/K 2 )/2. This the choice made in the 
chain calculation. 

(ii) By rescaling the variables as in section 3, one can convert the system into a homogeneous 
one with two defect lines. Repeating the transfer-matrix calculation, one finds again the spectrum 

CSP,flZ5p. 

(iii) The considerations here are not restricted to the critical point. However, the relation to 
the RDM via a conformal mapping is limited to criticality. 
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